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ABSTRACT 

We have examined the relationship between the optical depth of the 9.7 fim 
silicate absorption feature (rg.y) and the near-infrared color excess, E{J — Ks), 
in the Serpens, Taurus, IC 5146, Chameleon I, Barnard 59, and Barnard 68 
dense clouds/cores. Our data set, based largely on Spitzer IRS spectra, spans 
E{J — Ks)=0.3 to 10 mag (corresponding to visual extinction between about 2 
and 60 mag.). All lines of sight show the 9.7 fim silicate feature. Unlike in the 
diffuse ISM where a tight linear correlation between the 9.7 /im silicate feature 
optical depth and the extinction (Ay) is observed, we find that the silicate feature 
in dense clouds does not show a monotonic increase with extinction. Thus, in 
dense clouds, tqj is not a good measure of total dust column density. With few 
exceptions, the measured rg.y values fall well below the diffuse ISM correlation 
line for E{J — Ks)> 2 mag {Ay >12 mag). Grain growth via coagulation is a 
likely cause of this effect. 

Subject headings: ISM: clouds — ISM: dust, extinction — ISM: molecules: in- 
frared: ISM 
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1. Introduction 

Observationally, interstellar dust is known to consist mainly of amorphous silicates as 
evidenced by the 9.7 and 18.5 /im Si-0 stretching and bending vibrations of this material, and 
graphitic carbon, responsible for the strong 2175 A bump. Detailed models have been devel- 
oped that link the measured optical properties of these materials to the observed interstellar 
extinction using appropriate size distributions (e.g., Draine & Lee 1984; Desert, Boulanger, 
& Puget 1990; Zubko, Dwek & Arendt 2004). These models differ in their assumptions about 
the grain populations, whether silicates and graphite are physically separated or mixed in 
a composite grain with or without ices. Typically, these models conclude that both materi- 
als contribute about equal volumes (per H-atom) to the interstellar dust. The well-known 
Draine & Lee (1984) model illustrates that while silicates dominate the extinction beyond 
about 8 ixm through their strong resonances, in the visual through near-IR, graphite dom- 
inates the extinction. In contrast, some composite grain models argue that "big grains" 
made of silicates with a carbon-containing coating (Desert et al. 1990) dominate the extinc- 
tion in both the visible and IR. Zubko et al. (2004) consider a variety of dust components in 
their models including composite particles (containing silicate, organic refractory material, 
water-ice, and voids), polycyclic aromatic hydrocarbons (PAHs), bare silicate, graphite and 
amorphous carbon particles. The relative contributions of each of these components in the 
visible, near-IR and mid-IR varies and depends on the specific combination of components 
considered. 

In the diffuse ISM, observations spanning visual extinction {Ay) between 3 and 15 
mag have shown that the optical depth of 9.7 /im silicate absorption feature (tq 7) displays 
a tight linear correlation with Ay (Roche & Aitken 1984; Whittet 2003). This indicates 
that the silicate and graphite dust components are well-mixed and vary little in relative 
abundance. In dense clouds, however, there is observational evidence that this correlation 
suffers a significant break-down. The first indication that denser regions do not follow the 
diffuse ISM correlation came from the study of Whittet etal. (1988), who found that the 
correlation appeared to fail at high Ay. In that study, two (out of 5) lines of sight in the 
Taurus dark cloud showed rg.y values that were anomalously low with respect to their Ay. 

In this paper, we confirm and extend the failing of the rg.y versus extinction correlation 
for lines of sight in the IC 5146, Barnard 59, Barnard 68, Chameleon I, Taurus, and Serpens 
dense clouds. We show that a dramatic break-down of the diffuse ISM trend occurs when 
Ay exceeds ~ 12 mag. Our results are based primarily on Spitzer IRS spectra which we 
describe in f|2l In f|3] we discuss the relationship between the depth of the 9.7 /im silicate 
absorption feature, the visual extinction, and near-infrared color excess. Finally, in fQl we 
discuss the astrophysical implications of these observations in terms of grain properties in 
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dense clouds. 



2. Data Reduction and Analysis 

A summary of our dense cloud dataset including 2MASS IDs, spectral types, silicate 
optical depth at 9.7 //m, J — Kg, E{J — Ks), and Spitzer Program ID (PID) numbers, where 
relevant, are given in Table 1. The observations of the background sources behind IC 5146, 
Barnard 68 and Chameleon I were part of the Spitzer Cycle 1 program, PID 3320 (PI: Y.J. 
Pendleton). The sources were observed using the Infrared Spectrometer (IRS) modules Short- 
Low (SL) 2 (5.2-8.7 ;um), SLl (7.4-14.5 /xm), and Long Low (LL) 2 (14.0-21.3 /xm) providing 
resolving powers (A/AA) of 60-127 (SL) and 57-126 (LL2). These data were processed with 
the Spitzer pipeline version S12.0 to produce Basic Calibrated Data (BCD). Post-BCD 1- 
D spectra were created by carrying out background subtraction and custom extraction in 
IDL/SMART. The observations of the Barnard 59 sources were part of the Spitzer Cycle 2 
program, PID 20604 (PI: A. Boogert; T. L. Huard etal., in preparation). These data were 
processed with Spitzer pipeline version 15.3. The observations of the Serpens sources and 
the Taurus source, Elias 3, were obtained as part of the the "c2d" Legacy program (PI: 
N. Evans). The data reduction of these sources is described in Kncz ct al. (2005; and in 
preparation). Silicate optical depths for the Taurus sources Elias 9 and 13 are from Whittet 
et al. (1988). The silicate optical depth for Elias 16 is from Bowcy etal. (1998). We used 
the Spitzer IRS spectrum for Elias 3 as it has higher S/N than the spectrum presented by 
Whittet et al. (1988). Spectral types for all Taurus sources are from Whittet et al. (1988 and 
references therein). For the diffuse ISM, we plot data available in the hterature (see §3) and 
one additional measurement toward 2MASS 1743320-2847525 (t9.7=1.5, E{J - Ks)^4..21] J. 
Chiar, unpubhshed Spitzer data, PID 3616). 

To determine the optical depth of the 9.7 //m silicate absorption feature, we fitted a 
second or third degree polynomial across the combined IRS SL and LL2 (when available) 
spectra. The 5.2-7 and 13.5-15 jirci regions were chosen to represent the absorption-free 
continuum. In cases where LL2 data were available and the 18.5 //m silicate feature was 
very weak, the region between 20 and 21.3 /xm was also used as continuum. The optical 
depths were then calculated by T9.7 = — In Fgource /-^continuum, where F is the flux density in 
Jansky units. 

The J-Ks colors are from the Two- Micron All-Sky Survey (2MASS), except for 2MASS 
17112005-2727131, 17111538-2727144 (Barnard 59), and CK2 (Serpens). These sources were 
not detected in the 2MASS J band, and only an upper limit is given in the 2MASS catalog. 
For 2MASS 17112005-2727131, the source is also not detected in the H band, so we used 
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more sensitive ground-based H and Ks photometry from Roman-Ziifiiga et al. (2007) to 
estimate the J-band color. For all three sources, the J — Ks color is estimated by assuming 
median intrinsic colors appropriate for G0-M4 giants (Bessell & Brett 1988) and using the 
color-excess ratio E{J — H)/E{H — Ks) = 1.73 determined by Indebetouw etal. (2005) for 
lines of sight in the Galactic plane. 

All sources in Table 1 have 2MASS colors consistent with reddened background giant 
stars without intrinsic infrared excesses. In this way, we include only field stars located 
behind the dense cloud or core, and young stellar objects are excluded. In most cases, due 
to the intervening extinction, the field stars are not observable in the visible, and the visual 
extinction toward these stars cannot be directly measured. Thus, the near-IR differential 
extinction (color excess) is measured, and the visual extinction is inferred from that, based 
on an assumed extinction law. To eliminate the uncertainty that might be introduced by 
assuming a specific extinction law, we plot the near-IR differential extinction, E{J — Ks), 
in Figure 1. There is some uncertainty introduced for sources where the spectral types are 
not precisely known. This uncertainty is represented by the horizontal error bars in Fig. 1, 
and is not significant enough to change the trend of the data. 



3. The Relationship between E{J — Ks) and rg 7 

Figure 1 shows a plot of the measured optical depth of the 9.7 /im silicate absorption 
feature versus the near-IR differential extinction for the Taurus, IC 5146, Chameleon I, 
Serpens, Barnard 59 and Barnard 68 dense clouds. The solid line is the correlation line, 
Tg 7 = Ay/18.0, calculated based on 13 lines of sight in the diffuse ISM with Ay between 3 
and 29 mag (§2: Roche & Aitken 1984; Whittet 2003 and references therein), excluding the 
Galactic centecl. In order to plot this correlation on the rg.r versus E{J — Ks) plot, we 
converted Ay to near-IR differential extinction using the diffuse ISM extinction law given 
by Whittet (2003): Ay = 6.2x E{J — Ks). Comparing the loci of the dense cloud data with 
the diffuse ISM correlation line, we note the following. For E{J — Ks)< 2, the dense cloud 
sources follow the diffuse ISM correlation, showing a monotonic linear increase of rg.7 with 
E{J — Ks). Compared to the diffuse ISM trend, the dense clouds show a markedly weak 
silicate absorption feature relative to near-IR extinction for E{J — Ks)> 2 mag {Ay^ 12 
mag). With the exception of two lines of sight, rg 7 reaches a maximum value 0.85 for 



^Sources toward the Galactic center fall well above the correlation line due to their strong silicate absorp- 
tion features (t9.7= 3.2 — 4.0; Roche & Aitken 1985) relative to their Ay and color excess {E{J — K) ~ 5.0; 
Becklin et al. 1978). There, it is thought that the lack of carbon-rich dust in the central regions of the Galaxy 
reduces the extinction in the visual and near-IR relative to the extinction in the 10 /im region. 
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E{J — Ks) up to almost 10 mag (^y~ 60 mag). The outliers are two sources in the Serpens 
dark cloud. The Serpens source, SSTc2dl82852. 7+02824, is the only source in our dataset 
with E{J — Ks)> 2 whose silicate optical depth places it on the diffuse ISM correlation 
line. The highly extincted Serpens field star, CK2, stands out as having a very deep silicate 
feature, but it too falls well below the diffuse ISM correlation line. 

4. Discussion 

The available data show that in the diffuse ISM, the silicate extinction correlates linearly 
with E{J — Kg) (Whittct 2003). We present strong evidence that in dense clouds, this linear 
relation breaks down when E{J — Ks) exceeds 2 mag (equivalent Ay ~ 12 mag) (Figure 1). 
Thus, the data presented also show that r9.7 does not provide a good measure of total dust 
column in dense cloud environments. It is clear from Figure 1 that tq.t becomes weaker per 
unit E{J — Ks) as E{J — Ks) increases beyond 2 mag. The simplest explanation for this 
observed behavior is the effect of grain growth on the near-IR extinction. Grains of mean 
radius a produce extinction most efficiently when 27ra/A ~ 1 (Whittet 2003); for A ~ 2 //m, 
we have a ~ 0.3 /xm, which is a factor of about 2 larger than the grains responsible for visual 
extinction in the diffuse ISM. The effect is illustrated in Fig. 2, which plots opacities from 
Ossenkopf & Henning (1994) for the MRN size distribution (Mathis, Rumpl & Nordsicck 
1977, solid curve), compared with MRN subject to growth by coagulation (dashed curve) 
and coagulation + thin ice- mantle growth (dotted curve). Each curve is normalized relative 
to the opacity at 2.2 /im. The near-IR extinction is significantly enhanced by coagulation 
for a given silicate opacity. Ice mantle growth does not affect the strength of the silicate 
feature, but does affect the apparent shape due to blending with the H2O libration feature. 
This is expected since ice mantle formation adds only a negfigibly thin (< 200 A) layer that 
will not affect near-IR extinction much (Jura 1980; Draine 1985). Note that the shape of 
the 1-8 /xm continuum extinction is little affected by either growth mechanism (Fig.2). This 
is consistent with the observed invariance of near-IR extinction law in the Barnard 59 dense 
core over Ay between 6 and 60 mag. (Roman-Ziiniga etal. 2007). 

The various dust models proposed to explain interstellar extinction all assume that the 
extinction arises from a combination of silicate and carbonaceous grains. The main difference 
is in the assumed physical relationship between the silicate and carbonaceous components. 
In models that assume independent silicate and carbonaceous grain components, it is pri- 
marily the carbonaceous component that is responsible for the near-IR and visual extinction 
and the silicate grains that dominate the extinction in the 10 /xm region (e.g., MRN; Draine 
& Lee 1984; Kim, Martin & Hendry 1994). In contrast, models that assume multiple inde- 
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pendent components (small grains for the UV extinction, plus silicate cores covered with a 
carbonaceous mantle), the mantled silicate grains carry the extinction in the visible through 
the mid-lR (Desert et al. 1990; Li & Greenberg 1997). Similarly, in the models of Mathis and 
Whiffen (1989), the visible through mid-IR extinction is carried by a single grain component 
consisting of silicate and graphitic/amorphous carbon agglomerates. Finally, Zubko etal. 
(2004) considered a variety of dust components in their models including composite particles 
(containing silicate, organic refractory material, water-ice, and voids), polycyclic aromatic 
hydrocarbons (PAHs), bare silicate, graphite and amorphous carbon particles. They show 
that the contribution to the extinction in each wavelength regime is highly model dependent. 

Regardless of the details of the adopted dust model, grain coagulation in dense clouds 
appreciably affects the near-IR extinction. In dense clouds, significant grain growth is ex- 
pected to result from coagulation (Jura 1980). In the diffuse medium, grain processing is 
dominated by 50-100 km/s shocks. Because of the low strength of agglomerates, grain-grain 
collisions in even a weak interstellar shock will shatter agglomerates into their constituent 
"monomers" (Jones etal. 1996). Further processing, in the diffuse ISM, is then due to 
sputtering in shocks in the warm intercloud medium and reaccretion of gaseous silicon, mag- 
nesium, and iron atoms in an oxide mantle in diffuse clouds (Savage & Sembach 1996; Tielens 
1998). Thus, in diffuse clouds, these destructive processes prevent the grains from growing 
big enough to affect the near-IR extinction. Such a scenario could explain the good corre- 
lation of the near-IR/ visual extinction and the strength of the 9.7 /im feature in the diffuse 
ISM and the lack of a hnear correlation in dense clouds. Further laboratory studies on the 
coagulation behavior of mixed grain populations would be very helpful to settle this issue. 

Observationally, a number of questions remain: 

• Is the IR extinction law invariable in all dense clouds/cores? 

• At what extinction level or density docs the onset of rapid dust coagulation begin? 

• Does this onset vary from cloud to cloud? 

• Do ice mantles affect the coagulation efficiency? 

Further detailed studies that map the ncar-lR extinction and the 9.7 fim silicate absorption 
feature across the extent of individual dense clouds and cores are needed and will reveal 
valuable information about dust properties in these regions. 

This work is based [in part] on observations made with the Spitzer Space Telescope, 
which is operated by the Jet Propulsion Laboratory, California Institute of Technology under 
a contract with NASA. Support for this work was provided by NASA. The authors thank 
the referee, Viktor Zubko, for his helpful suggestions. 
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Table 1. Dense Cloud Lines of Sight 





Source 




opectral lype 


T9.7 


J — Ks 




opitzer rlU 


IC 5146:Qmdust 21-1 


21472204+4734410 


G0-M4 III 


0.51 


4.973 


4.16 


3320 


IC 5146:Quidust 21-2 


21463943+4733014 


G0-M4 III 


0.55 


2.733 


1.92 


3320 


IC 5146:Quidust 21-3 


21475842+4737164 


G0-M4 III 


0.50 


2.619 


1.81 


3320 


IC 5146;Quidust 21-4 


21450774+4731151 


G0-M4 III 


0.34 


2.185 


1.38 


3320 


IC 5146:Quidust 21-5 


21444787+4732574 


G0-M4 III 


0.33 


2.040 


1.23 


3320 


IC 5146:Quidust 21-6 


21461164+4734542 


G0-M4 III 


0.60 


4.337 


3.53 


3320 


IC 5146:Quidust 22-1 


2 1443293+4734569 


G0-M4 III 


0.55 


3.780 


2.97 


3320 


IC 5146:Quidust 22-3 


21473989+4735485 


G0-M4 III 


0.24 


1.197 


0.39 


3320 


IC 5146:Quidust 23-1 


21473509+4737164 


G0-M4 III 


0.50 


1.917 


1.11 


3320 


IC 5146:Quidust 23-2 


21472220+4738045 


G0-M4 III 


0.60 


2.474 


1.66 


3320 


Taurus ;EIias 3 


04232455+2500084 


K2 III 


0.40 


2.516 


1.76 


172 


Taurus:Elias 9 


04321153+2433380 


K2 III 


0.21'* 


2.109 


0.93 




Taurus:Elias 13 


04332592+2615334 


K2 III 


0.54" 


3.005 


2.26 




Taurus:EIias 16 


04393886+2611266 


Kl III 


0.75f 


5.443 


4.76 




Taurus :HD 29647 


04410804+2559340 


B8 III 


0.24<= 


0.597 


0.16 




Barnard 59 




17111538-2727144 


G0-M4 III 


0.85 


> 6.062 


5.58 


20604 


Barnard 59 




17112005-2727131 


G0-M4 III 


0.80 


> 4.409 


8.86 


20604 


Barnard 68 


Quidust 18-1 


17224500-2348532 


G0-M4 III 


0.2 


1.637 


0.83 


3320 


Barnard 68 


Quidust 19-1 


17224483-2349049 


G0-M4 III 


0.2 


1.720 


0.93 


3320 


Barnard 68 


Quidust 20-1 


17224407-2349167 


G0-M4 III 


0.45 


1.959 


1.15 


3320 


Barnard 68 


Velucores 1-1 


17223790-2348514 


G0-M4 III 


0.5 


2.938 


2.13 


3290 


Barnard 68 


Vclucorcs 1-2 


17224511-2348394 


G0-M4 III 


0.3 


1.440 


0.63 


3290 


Barnard 68 


Velucores 1-3 


17224027-2348555 


G0-M4 III 


0.4 


3.276 


2.47 


3290 


Barnard 68 


Velucores 1-4 


17224159-2350261 


G0-M4 III 


0.3 


4.044 


3.23 


3290 


Chameleon I: Quidust 2-1 


11024279-7802259 


G0-M4 III 


0.4 


1.701 


0.89 


3290 


Chameleon I: Quidust 2-2 


11055453-7735122 


G0-M4 III 


0.6 


3.375 


2.57 


3290 


Chameleon I: Quidust 3-1 


11054176-7748023 


G0-M4 III 


0.5 


2.361 


1.55 


3290 


Serpens: CK 2 


18300061+0115201 


mid K-carly M III 


1.8 


> 8.676 


9.45 


172 


Serpens: SVS76 Ser 9 


18294508+0118469 


G0-M4 III 


0.8 


4.326 


3.52 


172 


Serpens: SSTc2dl82852. 7-1-02824 


18285266+0028242 


G0-M4 III 


1.2 


4.136 


3.33 


172 



^Where spectral type range is given, it is based on placement in J ~ H versus H — Ks diagram, except for CK2. CK2 spectral 
type from Chiar et al. 1994. Taurus spectral types from Whittct ct al. 1988 and references therein. 

''2MASS photometry. Where lower limit is given, the source is not detected in the 2MASS J band due to its faintness. 

'^Computed using spectral types listed and associated intrinsic colors from Bessell & Brett 1988 

^Spitzer program ID for silicate spectra, where applicable. PID 3320, this work; PID 172, PI Evans, c2d Legacy; PID 3290, 
PI Langer. 

"Silicate measurement from Whittct ct al. 1988 

^Silicate measurement from Bowey, Adamson & Whittet 1998 



-10- 




Neor-IR Color Excess, E(J-Ks) (mog) 

Fig. 1. — Plot of the relationship between the optical depth of the 9.7 /im silicate absorption 
feature versus the near-IR color excess, E(J — Ks), for Taurus (triangles), IC 5146 (as- 
terisks), Serpens (diamonds). Chameleon I (squares), Barnard 68 (X's), Barnard 59 (open 
circles). Diffuse ISM lines of sight are represented by large filled circles (Whittet 2003, and 
references therein) and a small filled circle (Chiar, unpublished Spitzer data). The diagonal 
line represents the "diffuse ISM" correlation. Error bars indicate spectral type uncertainties 
in the calculation oi E{J — Ks). 
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Fig. 2. — Plot of dust opacities vs. wavelength. The solid curve is for the original MRN 
grain model applicable to the diffuse ISM. The dashed curve shows the effect of growth by 
coagulation after 10^ years at a particle density n = 10^ cm~^. The dotted curve shows the 
effect of thin ice-mantle growth in addition to coagulation. Each curve is normalized relative 
to the opacity at 2.2 //m. Data are from Table 1 of Ossenkopf & Henning (1994). 



